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Abstract 

The transformation of the partial fractional derivatives under spatial rota¬ 
tion in E? are derived for the Riemann-Liouville and Caputo definitions. 
These transformation properties link the observation of physical quantities, 
expressed through fractional derivatives, with respect to different coordi¬ 
nate systems (observers). It is the hope that such understanding could 
shed light on the physical interpretation of fractional derivatives. Also it 
is necessary to able to construct interaction terms that are invariant with 
respect to equivalent observers. 
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1 Introduction 


Fractional calculus deals with differentiation and integration to arbitrary 
real or complex orders. The idea is as old as the integer-order calculus. 
Extensive mathematical discussion of fractional calculus can be found in 
Refs. [l]-[5] and references therein. The techniques of fractional calcu¬ 
lus have been applied to wide range of fields, such as physics, engineer¬ 
ing, chemistry, biology, economics, control theory, signal image processing, 
groundwater problems, and many others. 

Physical applications of fractional calculus span a wide range of topics 
and problems (for a review see Refs. [6]-[9] and references therein). Gener¬ 
alizing fractional calculus to several variables, multidimensional space, and 
generalization of fractional vector calculus has been reported [10]-[13]. Also, 
progress has been reported on generalization of Lagrangian and Hamilto¬ 
nian systems [14,15]. 

Despite the significant progress in applying fractional calculus to a wide 
range of physical problems, there is still a lack of satisfactory geometric 
and physical interpretation of fractional calculus, in comparison with the 
simple interpretations of their integer-order counterparts (see Ref. [16] and 
references therein). Effort has been devoted to relate fractional calculus and 
fractal geometry [17,18]. A different approach to geometric interpretation 
of fractional calculus is based on the idea of the contact of a—th order [19]. 
However, a satisfactory interpretation is still missing [16]. 

A study of the symmetry of physical systems described by fractional cal¬ 
culus and investigating the transformation properties of fractional deriva¬ 
tive operators under specific groups can shed light on the nature of such 
operators. For example, studying the transformation of the fractional op- 
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erators under point transformation (or change of variables) can be found in 
Ref. [20, 21], even though the class of variables substitution preserving the 
form of the fractional derivative is found to be very narrow. Generalizing 
the definition of fractional derivative of a function with respect to another 
function can be used to study general transformation of point transforma¬ 
tion [22]. 

The question that we are trying to address in this work is basically 
the following: Given a quantity that is described by a fractional deriva¬ 
tive of some function in a specific coordinate system {x,y), how can this 
quantity be expressed in a rotated coordinate system {x',y')7 Thus, we 
try to relate the same fractional quantity in two equivalent coordinate sys¬ 
tems. In this work, we only consider the cartesian space coordinates in two 
dimensions, {x, y), and adhere to the notion of time as invariant in all co¬ 
ordinate systems (nonrelativistic). In other words, we are looking into the 
transformation of fractional derivative operators and their effect on physical 
quantities under the rotation group 50(2). We are not interested in con¬ 
structing the irreducible representations of 50(2), as done in the standard 
integer-order differentiation case [23], rather we focus on the orthonormal 
basis of the two dimesnional cartesian space of R^. Intuitively, quantities 
expressed through fractional derivative operators are expected to behave 
differently from scalar, vector, and tensor quantities. 

In this work we investigate the transformation properties of the frac¬ 
tional derivative and its action on an invariant scalar field, under space 
rotation in two dimensions, and using both the Riemann-Liouville and Ca- 
puto definitions. We compare between the Riemann-Liouville and Caputo 
definitions in relating the same physical quantity under rotation. We check 
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if there is a major difference in the transformation properties between the 
two definitions. In Ref. [24], the transformation properties of the Riemann- 
Liouville fractional derivative of a scalar field under infinitesimal transfor¬ 
mation of SO(2) are derived. The current work is more general and can be 
applied to any function and is not bounded to infinitesimal transformation. 
Also, we include the case of the Caputo definition of fractional derivative. 
Their special result agrees with our general results. In section 2, we give a 
brief introduction to fractional calculus to lay out our notation. In section 
3, we lay out the SO{2) transformation properties of fields and their deriva¬ 
tives expressed in Cartesian coordinates. In section 4, we apply our method 
to the Riemann-Liouville and Caputo fractional derivatives and layout their 
corresponding transformations under 50(2). Finally, in section 5 we give 
a brief discussion of our results. 

2 Fractional Calculus 

For mathematical properties of fractional derivatives and integrals one can 
consult Refs.[l]-[5] and the references therein. In this section we lay out the 
notation used in the rest of this work. We consider real analytic functions 
on and only discuss the Riemann-Liouville and Caputo definitions of 
the fractional derivative. Let f{x,y) to be a real analytic function in a 
specific domain in the Euclidian space R?] f : R? ^ R. The x-partial 
fractional derivative of order a (keeping y constant) and where a is the lower 
limit of X is written as aD'^f{x,y). Similarly for the y-partial fractional 
derivative of order a (keeping x constat) and where b is the lower limit of 
y is written as hDyf{x,y). Since the f{x,y) is analytic then the partial 
fractional derivatives are assumed to commute, i.e., Lt>“, bDy] = 0- 
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Definition 2.1. The Cauchy’s repeated integration formula of the nth- 
order integration of the function f{x,y) along x, keeping y constant, can 
he written as 


rx rXn-1 rxi 

aDf:'^f{x,y) = / dxn -1 dxn- 2 --- f{xo,y)dxo (1) 

J CL J Qi J CL 

^ 1 r f{u,y)du 

(n - 1)! Ja (x - 

A similar relation for the nth-order integration of the function f{x,y) 
along y, keeping x constant, 

ry rvn-i m 

bDy^f{x,y) =v dxn-i dyn- 2 ---J^ f{x,yo)dyo (2) 


1 


lb Jb 

fy f{x, u) du 


(n - 1)! Jb {y - 

Definition 2.2. The Riemann-Liouville fractional integration of order 
a < 0 and along x, keeping y constant, is defined as 


iDjfix,y) = 


1 


f{u,y)du 


r(-a) Ja (x - ’ 

and along y, keeping x constant, is 


bDyf{x,y) = 


1 


fy f{x,u)du 


( 3 ) 


( 4 ) 


T{-a)Jb ’ 

where r(.) is the Gamma function. 

We consider a = 6 = 0 and drop them from the notation henceforth. 
Definition 2.3. The Riemann-Liouville partial fractional derivatives of 
the order a > 0, where n — 1 < a < n and n € N, are defined as 


D^f{x,y) = Df:D^-^f{x,y) = 
D<Jf{x,y) = D’jD^-^f{x,y) = 


1 


f{u,y)du 


r(n - a) Jo (x - ’ 

1 fy f{x,u)du 
T{n - a) Jo (y - u) 


( 5 ) 

o-n+l ’ 


where to unify notation we used Df: = and Dlf = . 


— d" 
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The leibniz rule applied to the Riemann-Liouville partial fractional 


derivatives can be written as [l]-[5], 

r(a + 1) 


Dx {f{x,y)g{x,y)) = 


D<^-^f{x,y)Dlg{x,y), (7) 


^0 — k + l) 


D'i(f(x,y)g(x,y)) = f (8) 

Definition 2.4. The Caputo partial fractional derivatives of order a > 0, 
where n — 1 < a < n and n G N, are defined as 


C T^a 


Dxf{x,y) = 


1 


^Dyf{x,y) = 


_ r Dufiu,y)du 

r(n - a) Jo (x - ’ 

fy Df;^f{x, u) du 


(9) 

( 10 ) 


r(n - a) Jo [y - tt)“ ”'+^ 

The Riemann-Liouville and Caputo definitions of the fractional derivative 
are related [l]-[5] 


n—1 


X 


k — OL 


^Dff{x,y) Dff{x,y) g ^ ^ 

n—1 

^D^f{x,y) = D^f{x,y) \-y 


D\f{x,y)\x=o , ( 11 ) 


“ Dlf[x,y)\y=o . ( 12 ) 


As an example, for the function f{x,y) = x^y^, where /3,A > —1 we 


hnd 


Df (x%^) 

_ r(/3 + 1) ( B-a \\ 

- r{/3-a + l)y ^ ' 

(13) 




DfD'J 

r(/3 -h i)r(A +1) / a „ X 

T{/3-a + l)r{X-a + l)y ^ )’ 



The Caputo partial fractional derivatives give the same result for y^, 
except for f3,X G N < n, where in this case the Caputo partial fractional 
derivatives vanish, as expected from Eqs. illloD. 
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3 Transformation of Functions and Integer-Order 
Derivatives Under Spatial Rotation 

We consider the Euclidian space B? and work with orthonormal basis of 
the Cartesian coordinate system. The distinction between covariant and 
contravariant quantities disappears and we drop this notation in our discus¬ 
sion. As discussed before, we are not interested in discussing the irreducible 
representations of 50(2) [23], instead we keep working with the Cartesian 
coordinates (x, y) since we want to express the partial fractional derivatives 
in the same way. Consider the two Cartesian coordinate systems xy and 
x'y' sharing the same origin (0,0). The {x' ,y') coordinates are related to 
the {x , y) coordinates though a counter-clockwise rotation by an angle (j) 
about their identical axis of rotation z' = z. For (/> = 0 the two systems of 
coordinates are overlapping. In this work we consider the case of passive 
transformation of the coordinate system, thus we consider the observation 
of the same quantity from two different coordinate frames, without actually 
rotating physical quantities. The two coordinate systems are related by the 
transformation matrix 

Fcos^. ( 14 ) 

— sm (f) cos (p 

For the active transformation case the angle (f) would be replaced by (—</>)• 
The coordinates {x ', y') are related to the coordinates (x, y) through 



We present some of the properties of the functional transformations of 
the Lie group as applied to function f{x,y). The generator of the 50(2) 
rotation (angular momentum), around the z axis, L is written as, 

L = yDl - xD] . (16) 
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The group elements of the rotation group are written, in the exponential 
form, as where the exponential is defined as an infinite sum of differ¬ 
ential operators, explicitly 





(17) 


where (j) is the angle of rotation. Since Lx = y and Ly = —x, it is easy to 
derive the action of the element on the coordinates x and y, explicitly 


SL 

Qf rn - 


= COS 4>x + sin (f)y = X , 


y = — sin (px + cos (py = y'. 

Thus we retrieve the standard result in Eq. (1151) . 

Remark 3.1. In general, for n,m ^ Z~^ , 

x'^ y'^ = x”) y^) = (x” y^) . (18) 

This result emerges from proving the following two Lemmas. 

Lemma I: Given any two functions A{x,y) and B{x,y) that are analytic 
(of Class C°°) in their common domain, then 

n I 

= (19) 

A:=0 ' ' 

where n € Z~^. This is just the general Leibniz rule applied to the differen¬ 
tial operator L and the proof is similar using mathematical induction. The 
second important Lemma is 

Lemma II: Given any two functions A{x,y) and B{x,y) that are analytic 
(of Class C°°) in their common domain, then 

(^(x,7/).B(x,y)) = (e‘^^^(x,7/)) (e‘^^.B(x,y)) . (20) 
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Proof. We first write 


{A{x, y) B{x, 2/)) = £ —r {A{x, v) B{x, y)) 

r. 


( 21 ) 


n=0 

E (f) 71 TV ! 

1 un in — k)\ 

n=U fc=0 ^ ■' 






Noting that l/(n — A:)! vanishes for k > n we rewrite the above results as 


.<PL 


{A{x,y)B{x,y)) = J2J2 


(l‘a)( 


L^-’^B] . 


( 22 ) 


^ ^ k\{n - k)\ 

n=0k=0 ' ' 

Shifting the order of the sum and then shifting the index n —?> n + k, we 
write 


oo oo iri-\-k 

{A{x. y) B(x, v)) = E T, B) ■ (23) 

k=idn=—k 

Noting that 1/n! vanishes for n = —k, —k + 1,..., —1, we reach the final 
conclusion 


OO ik oo iji 

e*“-{A{x,y)B{x,y)) = '£^{r'‘A)'£^(r''m (24) 

1 ^ ' r\ 

fc=0 n=0 

= A{x,y)) B{x,y)) . 


Thus, Eq. (fTH]) follows. 

Based on Eq. (|18p we can discuss the transformation of functions un¬ 
der rotation. A function '^{x,y) that is analytic at (0,0), can be always 
expanded as 


^{x^y) = Y. 


n m 

I t ^ y ' 
ni ml 


In the x' y' coordinates frame, the transformed function ^'{x',y') can also 


be written as 


vI/'(x',y') = E 


i nn\ 


n\ m\ 


n,m 


x^y^- 
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If the constants a^m are independent of the coordinates, then we develop 
the following important result. 

Remark 3.2. A function ^(x,y) that is analytic in a specific domain, 
where its dependence on the rotation angle (p comes only through its argu¬ 
ments {x,y), is transformed as 

^'ix',y') = ^{x,y). (25) 

A scalar function ^{x,y) that is invariant under spatial rotation, must 
transform as 'I''(x',y') = 'I'(x,y), thus 

L'i>{x,y) = xDl'it - yDl'^ = 0 , (26) 

assuming the dependence on the rotation angle p occur only through the 
argument {x,y). 

Next we investigate the transformation of derivatives, one notes the 
following 

D]., x'^ = n = n e‘^^x^~^ = dI x^ (27) 

= Dl x^ = dI x'” . 

Thus we derive the important result 

01 = 6^^^016-^^ (28) 

and similarly for 

Ol, = 6"^^ Ol 6-'^^ . (29) 

The following commutation relations are easy to derive 

'oIl]=-oI, [oIl]=oI. 
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Transformation of higher derivatives can be derived, for example 

Dl = Dl dI dI . (30) 

In general we conclude 

Remark 3.3. The partial derivatives transform under spatial rotation as 
D^, = , D^, = , (31) 

where n,m £ N. 

As an illustration, one can show that 

Dl l] = -2DlDl , [d^l] = +2DlDl. (32) 

Thus, 

^1 + D^L = 0 , (33) 

indicating that the Laplace operator + Dy is a scalar operator, as ex¬ 
pected. 

An important remark is to check the transformation of the operator L 
itself under spatial rotation. 

D = L = L . (34) 

Thus the operator L is a scalar itself, a result that can be verified eas¬ 
ily. This is an important point as it allows for the definition of inverse 
transformation between the xy and x'y' coordinate frames. 
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4 Transformation of the Fractional Derivative Un¬ 
der Spatial Rotation 


In this section we generalize the results of the previous section to fractional 
calculus, in particular to the Riemann-Liouville and Caputo definitions of 
the fractional derivatives. We consider only real functions and real values 
of the fractional order a > 0. This restriction limits the allowed domain 
space {x, y). This is clear if we consider the transformation of and , 
where P £ R. 

^<f>L ^/3 _ (j)x + sin (j) y)^ = x'^, 
g<t>L yP _ 0 2 ; + cos (j) y)^ = y'^ ■ 

Therefore, we require x, y > 0 and x', y' > 0, which is equivalent to only 
consider the allowed domain: [(x,y) : x,y > 0 fl yjx > tan(/)]. Despite the 
restriction on the domain, one can still study the transformation of the frac¬ 
tional derivatives and deduce their properties. An extension to the whole 
domain can be pursued in the complex plane once a clear understanding of 
such transformation is acquired in the restricted domain. 

Based on the results of the previous section, we have 
Remark 4.1. In the domain [x,y > 0 fl yjx > tan(/)], the quantity x^y'*' 
transforms as 


x'^ y'^ = (e‘^^ x^ 




(35) 


where I3,X £ R. Next we look at the transformation of the partial frac¬ 
tional derivatives. Using the Riemann-Liouville and Caputo definition of 
the fractional derivative we note the following 


D“,x'^ = 


r(/3 + 1) ,y-a 
r(/3 — a -|- 1 ) 


= x^ = x^ = D" x'^ 


T{(3-a + l) 


(36) 
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Remark 4.2. The partial fractional derivatives according to the Riemann- 
Liouiville and Caputo definitions transform as 

(37) 

Dp = Dp . 

Considering an infinitesimal angle transformation, (p «1, and keeping 
only linear terms of (p, we write 

DP^DP + P[L,DP] , (38) 

Dp^Dp + p[L,Dp . 

Remark 4.3. Given an analytic function f{x,y), its partial fractional 
derivatives transform as 

DPf\xPy') = e^^DPf{x,y), (39) 

Dpf'{x',y') = e^^D^f{x,y). 

Considering the infinitesimal angle approximation, p « 1, and keeping 
only linear terms of p, we write 


DP fix', y') . 

- Dpf{x, y) + pLDpf{x, y ), 

(40) 

Dpf'{x',y') 

~ Dyf{x, y) + pLDpf{x, y). 

(41) 


Thus we derived the transformation of the partial fractional derivatives ac¬ 
cording to the Riemann-Liouville and Caputo definitions. Next we consider 
the transformation of the partial fractional derivatives as applied to a scalar 
field T(x,y) and using both Riemann-Liouville and Caputo definitions. 
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4.1 Transformation of the Riemann-Lionville Partial Frac¬ 
tional Derivatives of a Scalar Field 

Consider a scalar function ^{x,y) that is invariant under spatial rotation, 
i.e., L'I'(x,y) = 0, then according to Eq. ()10]l 


y') « y) + y ), 


(42) 


where 


LZ)“4/'(x', y') = yDlD^-^{x, y) - y). (43) 

Since y and Dy commute with 1)“ we can write the above equation as 


LD^^'ix', y') = (2/4/) - . (44) 

Rewriting the above equation using commutation relations, we find 
LZ?“4/'(x', y') = [dI, (2/4/) - [x, (z)^4/) + (yDl - xD]) 4/ .(45) 


The last term vanishes since L4/ = 0. Thus we arrive at 
LZ)“4/(x,2/) = [dId<^] (2/4/) - [x,Z)“] (z)i4/) 


(46) 


Note that from the properties of the fractional derivative, = D™'*'", 

while [DI.,D^] ^ 0 except for a G . We use the known relation [l]-[5] 

(4/(0,2/)) 


D<iDl^{x,y) = D<^,+^^{x,y)- 
Thus, we conclude 

(2/4'(x,2/)) = 4 




n-a) 

2/4'(0,2/)x“““^ 


(47) 


(48) 


r(-a) 

Using Leibniz rule in Eq. ([7]), substituting /(x) = 4/(x,2/) and g{x) = x, 
one can show that 


[x,D'^]Dl^!{x,y) =-aD'^ ^T)^4'(x,2/) 


(49) 
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Thus we conclude that 




y'^{o,y)x " ^ 


^ y r(-a) 


(50) 


A similar relation holds for y 






—a—1 


r(-«) 


(51) 


Note that for a G Z'*', the last term in the above result vanishes. Also for 
0 < a < 1, the term Dy~^ represents a fractional integral 
Remark 4.1.1 For a = 1 we recover the standard transformation of the 
gradient operator, 


Dl,^'{x', y') = Dl^{x, y) + (t^ Zl^T(x, y) , (52) 

Dl,^'{x',y) = Dl'^{x,y) - (/)Il^T(x, y). 


Remark 4.1.2 For a = 2 we recover the known scalar invariance of the 
Laplacian operator 

(dI + Dl) T'(x', y') = (dI + Dl) T(x, y ). (53) 

4.2 Transformation of the Capnto Partial Fractional Deriva¬ 
tives of a Scalar Field 


To derive the transformation of the Caputo partial fractional derivatives 
we consider 0 < a < 1, and use the relation between the Riemann-Liouville 
and Caputo definitions, as expressed in Eqs. (HHIll). Also we use the 
results of Eqs. One can show that 


x, 


Cd^ 




x^-^Dl^{Q,y) 
F(1 - a) 


(54) 


'T' 1 T 


iy^)=\DlD^] (yT)- 


X “ ^y'I'(0, y) 

rra 


(55) 
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Therefore, the Caputo partial fractional derivative transform, under in¬ 
finitesimal angle 4> « 1, as 


^ ^D^-i>{x,y) + cpLD^-^Dl^+ '' • ( 56 ) 


1 “ n,'\ ~ C n«lTt/™ J. / 1 7~ll,Tt I y 




r(l - a) 


• (57) 


Note that a — 1 < 0 and thus ^ and Dy ^ represent fractional integrals. 


5 Discussion and Conclusions 


In section 4, we derived the transformation of the Riemann-Liouville and 
Caputo partial fractional derivatives under spatial rotation. We also de¬ 
rived the transformation of the Riemann-Liouville and Caputo partial frac¬ 
tional derivatives of a scalar function. We observe no fundamental dif¬ 
ference between the transformations of the Riemann-Liouville and Caputo 
partial fractional derivatives. The fact that the Caputo fractional deriva¬ 
tive of a constant vanishes does not lead to a significant difference in its 
transformation. The difference we observe is that in the Riemann-Liouville 
transformation, there is a term that depends on the value of the function 
at the boundary. While, for the Caputo transformation, a similar term de¬ 
pends on the normal derivative of the function at the boundary. We observe 
that the Riemann-Liouville and Caputo partial fractional derivatives of a 
scalar function do not transform as the components of a gradient vector, as 
in the case of the integer-order, a = 1, partial derivatives. Thus they are 
classified to have a new form of transformation, different from the scalar, 
vector, and tensor quantities. They represent new entities that deserve to 
be studied thoroughly in future work. 
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Next we provide few examples of scalar functions. Consider the con¬ 
stant scalar function f{x,y) = 1, According to Eqs. (|^[^T|) . the Riemann- 
Liouville transformations of the fractional derivatives are, for infinitesimal 
(j) « 1 , 






y X 


— Q.— 1 


T{-a) ^ 

^ — Ci — 1 

xy 


(58) 


r(-a) ) • 

One can easily check that the result is as expected, for example 


Dm = 


x' " _ {x + <pyy 

r(-a) r(-a) 




yx 


-0.-1 


n-a) ’ 


(59) 


as required. It is easy to show that the following quantities are invariant. 


x“Z)“(l), y“D“(l) (60) 

y-“ZI“(l) - X-“Z)“(1) , x2+“ZI“(l) + . 


One can use {dx)°‘ and {dy)°‘ instead of x" and ?/“ in the above invariant 
forms. The Caputo partial fractional derivatives of a constant field vanish. 

Consider the scalar function 'I'(x, y) = x'^ + The Riemann-Liouville 
and Caputo x-partial fractional derivatives are, respectively 


T)"T(x,y) 

^Z)“T(x,y) 


2-a I 2 X 

-X +y 


r(3-a)'^ ' " r(i - a) ■ 

2 


r(3 — a) 


X 


2-a 


(61) 


Again it is easy to check that the transformations of the partial fractional 
derivatives, for infinitesimal angle (j) « 1, agree with Eqs. (|50ll56p for 
both the Riemann-Liouville and Caputo fractional derivatives. A similar 
transformation holds for and ^Dy with x and y interchanged and cj) —>■ 
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—([). Once can easily check that the quantities 

y) + y), x'^y^D'^ y ), (62) 

x“ ^D'^^ix, y) + y). (63) 

are invariant under spatial rotation. Again, one can use (dx)^ and (dy)“ 
instead of x" and in the above invariant forms. 

For the function 'l'(x, y) = (x^+y^)^, one can check the transformations 
correct. However, the quantities in Eq. (1621) are not invariant separately. 
Nevertheless, it is possible to find a linear combination of the two quantities 
to form an invariant quantity, namely, 

x“Z)“M/(x, y) + y“D“'k(x, y) + A [x^y^D^ Z)“'I/(x, y)) . (64) 

where A is some constant to make the combination invariant. Again, one 
can use (dx)“ and ((iy)“ instead of x“ and y°‘ in the above invariant forms. 

In Ref. [24], the authors derived the inhnitesimal transformation of the 
partial fractional derivatives of a scalar held, by using the function’s Taylor 
expansion and then applying the Riemann-Liouville fractional derivative to 
each term of the expansion. Their result agrees with our general result. 
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